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We study the anomalously large specific heat jump and its systematic change with pressure in 
CeCoIns superconductor. Starting with the general free energy functional of the superconductor 
for a coupled electron boson system, we derived the analytic result of the specific heat jump of the 
strong coupling superconductivity occurring in the coupled electron boson system. Then using the 
two component spin- fermion model we calculate the specific heat coefficient C(T)/T both for the 
normal and superconducting states and show a good agreement with the experiment of CeCoIns. 
Our result also clearly demonstrated that the specific heat coefficient C(T)/T of a coupled electron 
boson system can be freely interpreted as a renormalization either of the electronic or of the bosonic 
degrees of freedom. 

PACS numbers: 74.20, 74.20-z,74.50 



The superconductivity (SC) in CcMIn 5 (M=Co, Rh, 
Ir) heavy fermion (HF) series compounds 1,2 provides 
valuable information about the superconductivity nearby 
a possible quantum critical point (QCP) due to the excel- 
lent tunability with pressure and magnetic fields 3 . Be- 
sides a complicated phase diagram due to a compete- 
tion/interplay between the magnetism and superconduc- 
tivity, the superconducting properties as well as the nor- 
mal state ones around the QCP exhibit various anoma- 
lous behaviors such as non- Fermi liquid (NFL) power 
law of temperature in the resistivity, pseudogap behavior 
above the superconducting critical temperature, the con- 
tinuous evolving from the second order to the first order 
phase transition of the superconductivity with magnetic 
fields, etc 3 . 

In a recent paper 4 we proposed the two-component 
spin-fermion model as a generic pairing mechanism of 
this series of compounds CcMIns. In this work, it is 
shown that some of the standard features of the d-wave 
superconductivity mediated by magnetic fluctuations be- 
comes strongly modified in the vicinity of the magnetic 
critical point. Indeed it is very natural to expect that 
such a strongly coupled fermion-boson superconductiv- 
ity shows various anomalies deviating from the standard 
weak or strong coupling theory of superconductivity 5 . 

In this paper we study anomalous property of the spe- 
cific heat coefficient C(T)/T in CeCoIns. Sparn et al. 6 
measured C(T)/T of CeCoIns with varying pressure and 
also with applying magnetic fields. The key findings are 
that (1) suppressing the superconductivity with 8 Tesla 
fields C(T)/T show InT divergence down to 0.4 K; (2) 
without the magnetic fields the same sample becomes 
superconducting at 2.3 K with the specific heat jump 
ratio AC(T)/C(T) ~ 5, which is a huge jump consider- 
ing C(T C )/T C ~ 400mJ/K 2 mol; (3) with applying pres- 
sure T c gradually increases to 2.8 K at the pressure of 15 
kbar but the jump ratio AC(T)/C(T) drops drastically. 
We think this peculiar behavior is another manifestation 
of the critical magnetic fluctuations mediated supercon- 
ductivity. The huge specific heat jump, indicating the 



steep drop of the entropy below T c , is actually consistent 
with the spin-fermion model if we consider that the crit- 
ical magnetic fluctuations in this model is generated by 
the fermionic particle- hole excitations, which should be 
gapped below T c . This means that below T c not only the 
fermions suffer the entropy drop by developing pairing 
correlation but also the spin fluctuations should suffer the 
entropy drop, which add up to make a larger specific heat 
jump. As will be shown below, however, this phenomena 
is actually very general so that it should occur even in 
the electron-phonon mediated superconductors. The rea- 
son why it is phenomenal in CeCoIns is simply because 
the coupling of fermion-boson is strong and in particu- 
lar the resulting bosonic mode is near critical. There is 
a recent work 7 addressing the same problem as in this 
paper, but proposing a different mechanism for the en- 
hanced AC(T)/C(T). 

We start with the general free energy functional of the 
two component spin-fermion model, which is nothing but 
the same form as the one for electron-phonon model 8-10 . 
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where P = (p, iu n ), q = (q, iO„); u n = nT(2n + 1), fl n = 
7rT(2n), and G(P) = (iw n + e p + £i(-P))/<£(P), P(P) = 
-£ 2 (-P)/>(P), 4>(P) = (»w„ - Cp - £i(P))(^„ + e p + 
Ei(— P)) - |£ 2 (P)| 2 , respectively. Finally, the dressed 
spin propagator is written as D~ 1 (q) = D^ ) 1 {q) — Tr{q). 
flf er , £l sp i ni and O int are just apparent separation of the 
total free energy but the literal meaning of them is not 
necessarily correct as will be clear later. The notations 
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are standard 10 , and the one important difference is that 
Do(q) is not explicitly defined in the spin-fermion model 
in contrast to the electron-phonon case, but only the total 
spin propagator is defined phcnomcnologically as follows. 
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where I{T) defines the distance from the magnetic quan- 
tum critical point (QCP) and A and T are the coeffi- 
cients determined by the Fermi liquid parameters of the 
fermion sector 11 . The effective spin correlation length 
and the relaxational energy scale are £ 2 = A/I(T) and 
uj s f = r • I(T), respectively, and the momenta q and Q 
are assumed to be in the two dimensions 12 . We take V as 
the unit energy scale in our numerical calculations. The 
bare spin propagator D is assumed to be trivial regard- 
ing the critical low energy behavior and taken to be a 
constant. As usual, the above free energy functional is 
constructed in such way that it is stationary with respect 
to the variations of Ei^ and ir a la Luttinger-Ward 8 with 
the following definitions. This property of the functional 
is very important and utilized critically in our paper. 
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This is our key result in this paper and we can read 
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the specific heat jump ratio as 

We see from this formula that the jump ratio is modified 
with the normal self-energy renormalization factor 
which is i in the BCS limit with Z = 1 and can have a 
maximum value 1 in the limit of Z(T) — > oo. Apparently 
t^Pj- — > 1 is the proper limit for CeCoIn 5 and some other 
heavy fermion compounds exhibiting a large C(T)/T at 
low temperature. We summarize the consequences of this 
formula: (1) The jump ratio at most can be doubled to 
the BCS ratio assuming Q^p- is the same as the BCS limit, 
i.e. 9.42T C ; (2) Apparently the steeper slope of A(T) can 
contribute an extra enhancement to the jump ratio. And 
since our formula is an exact result, this formula can be 
utilized to estimate the temperature slope of A(T) at T c 
from the experimental in the limit of —> f; 

(3) Finally, the interpretation of the large C(T)/T and 
its enhanced jump ratio from the heavily renormalized 
fermionic quasiparticle due to the scattering from the 
magnetic fluctuations seems to be legitimate, if we can 
neglect the contribution from r2f, oso „ 13 . 

Now in order to calculate Ctotai(T) above and below 



T c , we need to calculate j with a given Z?(q)(")'( s ' . 
(5) But we found a more convenient way to proceed as fol- 
lows. Adding flint term to f2 sp m exactly cancels the cum- 



bersome term n(q)D(q) as 



n spm + n mt = (-T/v^lH-D-'iq)]. 
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Now following Ref[ 10 ], neglecting the momentum de- 
pendence of £^2 and using the definitions, Ei(-P) = 
iu) n (l — Z(u> n )), and A(w n ) = T,2(P)/Z(u! n ), we obtain 
after the momentum integration the following. 



n fer + fl int = -T.N(0)-n 
A 2 
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Note that the above equation didn't include the £l S pin 
term yet 13 . But due to the stationary property of the to- 
tal free energy functional as mentioned above, the specific 
heat jump across the superconducting transition with the 
changes of E^™' ) — > and — > is solely deter- 
mined by Eq.(7) with Z(u) n ) = Z^(ui n ) at T = T c . 

The Eq(7) can be written after the Matsubara fre- 
quency summation using the contour integration as 
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where uj 2 = e 2 + A 2 and the specific heat is calculated at 

T = T r as 



Furthermore, f£/ er term alone, after the same manip- 
ulations used for Eq.(7) and Eq.(8), becomes 



fl fer = -T ■ N(0) ■ ttY^PVA 2 -(iu n ) 2 - 2 
= 7V(0)^ de 2lj ■ tanh(^). 
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A different organization greatly simplifies £1 f er so that 
all the renormalization effect due to the interaction 
miraculously drops and furthermore the explicit A(T) de- 
pendent term all cancels out in the final result of Eq.(ll). 
This amazing result is actually the consequence of the 
Eliashberg free energy functional built with the station- 
ary property with respect to the variations of the self- 
energies. Considering that the Eliashberg free energy 
functional is not the exact free energy functional but a 
Consistent approximate functional equivalent to the self- 
consistent Born approximation for the fermionic part and 
the random phase approximation (RPA) for the bosonic 
part, the exactness of Eq.(9) and Eq.(ll) has a limited 
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meaning. Nevertheless, besides the accuracy of this func- 
tional the insightful reason for this result is the equivalent 
relation of the following. 
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From the technical viewpoint the above relation is triv- 
ial if we consider any linked-cluster diagram obtained 
from the free energy expansion with Hi nt s . However the 
physical interpretation is rather revealing. Namely, the 
interacting free energy Afi = Cltotal ~ Qbare and its de- 
rived entropy AS all comes from the expansion of Hi nt . 
Then depending on how to view or organize the same 
diagrams, we can view all interaction effect as the renor- 
malization of cither the fermion sector or the bosonic 
sector 15 . While theoretically we have this freedom, in 
experiments the experimental probe itself determines 
how to measure the contribution of Hint through either 
fermionic or bosonic degree of freedom. However the 
thermodynamic measurement like the specific heat mea- 
sures simultaneously the both degrees of freedom and it is 
hard to know how much contribution comes from which 
degrees of freedom. Our analytic result clearly demon- 
strated that the specific heat coefficient C(T) /T of a cou- 
pled fermion-boson system can be freely interpreted as a 
renormalization either of the electronic or of the bosonic 
degree of freedom, but not for both. This point should 
have a far-reaching consequence to the heavy fermion ex- 
periments and its related interpretations. This important 
aspect will be dealt in the separate publication. 

Now let us return to the calculation for the total 
C(T)tot- From Eq.(lO) we can calculate the entropy 
S(T) sp i n+ i nt with the phenomenological D(q) given in 
Eq.(3) as follows. 
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In real calculation, we need to fix parameters. Our 
unit energy scale is T = 1, and we take A w = 3, 
A(q' nax - Q) 2 = = 1.0, and a = 1.0. For the 

fermionic part f2/ er , it is just the same form as the non- 
interacting fermion free energy except lu = \/ t'f, + A 2 (T), 
and therefore it develops some structure below T c but no 
discontinuity at T c . Finally to fix the relative magnitude 
of rt sp i n+ i nt to Slfer wc need to determine the spin den- 
sity of states N s defined as N s f^" = y J^*" * . This 
can be determined using Eq.(12) at normal state, which 
gives the following relation 



i(T) + uj™ ax /r + \n n \/T 
/(T) + |Q n |/r J 



We assume n n = — 



|j]„| 



in the derivation, but the above 



relation is not useful since both side of the equation di- 
verges unless we introduce high frequency cut-off for Q n 
and io n - Since the above relation should hold for a vari- 
ation of <57r n , we derive another equivalent relation at 
T = as 
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A recent paper by Ref.[ 16 ] also use the similar pro- 
cedure to fix the relation between N(Q) and N s . Our 
numerical result gives N s = 1.5 ^ A^(0). Now we are 
ready to calculate C to t(T) using Eq.(ll) and Eq.(13). For 
the calculation below T c we need to know the tempera- 
ture dependence of the gap function A(T) and 7r s (A(T)). 
For the gap function wc use the generalized BCS form 
A(T) = A tanh(a y/T c /T - 1), where a = 1.764 and 
Aq = 1.74T C for BCS limit, but for strong coupling su- 
perconductor these values can deviate largely from the 
BCS limit and wc take them as parameters. Now much 
unclear part is the behavior of 7r s (T) below T c . Quali- 
tatively it should be cut off for fl < 2Aq at T = and 
smoothly recover to the normal state form approaching 
T — > T c . The leading expansion of 7r s (T) in A 2 (T) gives 
7r s = 7r„ cxp(-A 2 (T)/[w 2 + T 2 ]); but for a larger value 
of A(T) for T — > this form should not be very correct. 

The numerical results are shown in Fig.(l) and Fig. (2). 
In Fig.(l) the specific heat coefficient C(T)/T calculated 
from Eq.(13) are shown in unit of D °S N(0). For our pa- 
rameter choice producing Z >> 1 (For this we need the 
critical spin fluctuations (Jo << 1) as well as the effective 
coupling of 0(1)), the contribution of C{T) /T from 
tlfer Eq.(ll) is negligibly small. Therefore we didn't 
add it to C{T)/T in Fig(l) but show it separately in 
Fig. (2) demonstrating its qualitative features. In Fig.(l) 
the solid line is for the normal state C{T)/T with Iq = 
(the spin fluctuations is at QCP), it indeed displays InT 
divergence with decreasing temperature. With SC tran- 
sition at T=0.2 (open square symbol) it shows the jump 
ratio AC(T)/C(T) ~ 5 with the choice of parameters 17 . 
Then increasing Iq, the InT divergence is quickly sup- 
pressed and at the same time the specific heat jump ra- 
tio drops with increasing T c (T c is increased by hand and 
the temperature slope of A(T) is reduced accordingly as 

. /' { = tf 2 -)- With this rather crude phcnomenologi- 

A (T c2 ) Jci ' 1 ° 

cal choice of parameters our numerical results reproduce 
the experimental features of C(T)/T in CcCoIns quali- 
tatively as well as quantitatively. 
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In summary, starting with the general free energy 
functional for the coupled fermion-boson system, we de- 
rived an analytic formula of AC(T)/C(T) for the gen- 
eral strong coupling superconductor. Then we calculate 
C(T)/T for the spin-fermion model and show that the 
salient features of C(T)/T of CeCoIn5 for both normal 
and superconducting states are successfully explained in- 
cluding its anomalous jump ratio AC(T)/C(T) and the 
progressive reduction of it with increasing T c . Also with 
different organization of the Eliashbcrg free energy func- 
tional, which lumps all the interaction effect cither into 
the bosonic or into the fcrmionic degrees of freedom, we 
clearly demonstrated that the effect of the interaction in 
the total free energy or entropy of the coupled fermion- 
boson system can be freely viewed as the renormalization 
of either the fermion sector or the boson sector. Finally, 
our results strongly support the idea that the two di- 
mensional critical magnetic fluctuations plays an essen- 
tial role in CeMIns HF compounds in producing a large 
and strongly temperature dependent C(T)/T and the SC 
pairing itself. 
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